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SUMMARY 


Two methods are presented for dynamically scaling store models for 
wind-tunnel store -dropping studies. For each method the model and proto- 
type Mach numbers are assumed to be necessarily equal and the Reynolds 
number effects are assumed to be negligible. The light-model method gives 
exact simulation of the store motions except that the vertical displacements 
are deficient. This deficiency is reduced as the vertical ejection veloc- 
ity is increased or it can be eliminated by accelerating the parent model 
upward at the instant of store separation. The heavy -model method, in 
which the parent model is stationary, gives complete simulation of the 
store motion except that the short-period longitudinal oscillation is too 
poorly damped; this defect is of no serious consequence during the first 
phase of a drop because the time during which the store is critically 
close to the parent model is generally small compared with the period of 
the longitudinal oscillation. The heavy -model method is generally recom- 
mended for store-dropping studies; however, it is often impossible to 
make the models suff iciently heavy and with the proper moment of inertia. 

If such is the case, the light-model method is required. 

A brief description of the method of conducting store-dropping tests 
in the preflight jet of the Langley Pilotless Aircraft Research Station 
at Wallops Island, Va., is given. • 


INTRODUCTION 


Investigators have, for some time, been using the free-model tech- 
nique in wind tunnels to study the motions of stores and missiles during 
and following release from a parent airplane. In this method, a model of 
all or part of the parent airplane is placed in the tunnel test section 
and the store or missile model is dropped or propelled into the airstream. 
The unrestrained dynamic behavior of the store or missile model in the 
presence of the interference field is recorded by means of high-speed 
photography . 
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In making tests in the aforementioned manner, it is necessary, of 
course, to scale the models in a manner which will determine the relations 
between the translational and rotational motions of the model and the 
prototype. The general subject of scaling aerodynamic models for free- 
flight testing has been thoroughly discussed in references 1 and 2. In 
the present report, the scaling relations have been derived in a manner 
which permits an insight into the specific problems of free-model drop 
testing in wind tunnels. Some of the deficiencies of the test method 
are noted and methods by which the initial conditions may be adjusted 
to circumvent these deficiencies are presented. In addition, a brief 
description of the free-flight drop-model technique as employed in the 
preflight jet of the Langley Pilotless Aircraft Research Station at 
Wallops Island, Va., is given in the appendix. 


SYMBOLS 


x 

Cx 


7 

P 

M 



linear horizontal displacement 

aerodynamic force coefficient in horizontal plane 
ratio of specific heats 
static pressure 
Mach number 


S 

L 

N 

d 

g 

% 

V 



reference area 

characteristic length 

number of characteristic lengths 

average model weight density 

acceleration due to gravity 


volume 
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time 

linear vertical displacement 


K z = 

Cz 

tx 

0 

C0 

rr _ k 

k R - 7 


c z i 


K V 


aerodynamic force coefficient in vertical plane 

time required to move a specified number of characteristic 
lengths in horizontal plane 

pitch angle 

total aerodynamic moment coefficient 


KttK, 


fc V R 


k radius of gyration 

r c e I 

K 0 ~ — 

a z,tx 
U 
P 

C % 

t l/2 
W 

q. 


( Cm e + Cm i) 


angle of attack contributed by z at t 

stream velocity 

period of oscillation 

static pitching -moment derivative 

time to damp oscillation 
weight 

dynamic pressure 

total damping derivative 


= 


J a 


lift-curve slope 
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T 


stream static temperature 


*ax 


specified horizontal load factor 


Kaz 


specified vertical load factor 


Subscripts : 


0 


condition at t x = 0 


M 


model 


P 


prototype 


Dots over quantities indicate derivatives with respect to time. 
The coordinate system is illustrated in figure 1. 


In making free-model tests, it is necessary to scale the model 
dimensions, weight, and inertia in a manner which will determine the 
relations between the translational and the rotational motions of the 
model and the prototype. The primary requirement is that the trajectories 
of the model and prototype be as nearly geometrically similar as possible. 
The method of analysis consists of writing expressions for the transla- 
tional and rotational motions and for the period and damping of the short- 
period longitudinal oscillation. Examination of these equations indicates 
those parameters required for the particular simulation desired. Two 
basic methods of simulation are described; in each, it is assumed that the 
model and prototype Mach numbers are equal and that Reynolds number effects 
are negligible. In the first method, the aerodynamically produced trans- 
lational accelerations vary inversely with the characteristic length; this 
method is best suited for dynamic stability studies and for store-dropping 
studies where a large vertical ejection velocity is used. In the second 
method, the aerodynamically produced accelerations are independent of 
characteristic length just as are the gravity-produced accelerations; this 
method is required for free -dropping store studies. The first and second 
simulation methods will be called the light-model method and the heavy- 
model method, respectively, for the purposes of the present report. 

The present method of analysis may appear more tedious than the 
classical dimensional analysis given, for example, in reference 3- Both 
methods yield the same dimensional relations between the model and the 


MODEL-SCALING EQUATIONS 
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prototype. However, the present method shows how some of the deficiencies 
inherent in the model tests may be circumvented by adjusting the initial 
conditions . 


Light-Model Method 

The equations for linear motions in the horizontal plane are 


C x | pM 2 K g L 2 g K^pg 
dKylr'* ” dL 


X = X 0 + 


KxPg 

dL 


( 2 ) 


X 


x 0 t + 


+2 

2dL % 


(3) 


where xq is the horizontal velocity at t = 0 and 
quantities which are assumed to be independent of L. 


K x contains those 
The quantity 


dL 


is the ratio of aerodynamic force to weight (load factor); if ^ 


is constant, the load factor and, consequently, x 
tional to L. 


are inversely propor- 


An expression for the time t x required to move a specified number 
of characteristic lengths NL in the horizontal plane is required. From 
equation (3) with xq = 0 and x = NL, the following is true: 



00 


Therefore, t x is proportional to L if 2. q s constant. Substituting 
t x (from eq. (4)) for t in equations (2) and (3) shows that x^ x and 
x tx (conditions when t = t x ) are independent of and proportional to L, 
respectively. 
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The equations for linear motions in the vertical plane are 



( 5 ) 


i 


z 0 + 



( 6 ) 


Z = ZQt + 



K zP\ gt 2 
dL J 2 


(7) 


where 



is the vertical load factor. 


Substituting t x (from eq. (4)) for t in equations ( 6 ) and ( 7 ) 
shows that z and z approach independence of and proportionality to 
L, respectively, only as the vertical load factor or Zq becomes large. 

Consequently, the vertical and horizontal motions do not scale identically 
for the conditions so far established. The vertical launching velocity 
may, however, be adjusted so that the ratio of vertical to horizontal 
displacement at a specified value of t x will be equal for model and 
prototype. Substituting t x (from eq. (4)) for t in equation ( 7 ) and 
dividing by x = NL yields 




dL 



(8) 


It is desired that this ratio be equal for model and prototype when 
t = t x as defined by equation (4). Writing equation (8) for model and 

prototype gives 



(9) 


where subscript M denotes model nnd subscript P denotes prototype. 
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proto - 

factor 
0.5 foot 



= ZQ^p + 12.1 ft/sec 


This value of z will make ( — ] equal for model and 
°' M \ x /t v 

KyP 

type. For a typical case where is such as to give a load 

of 20 for a 20 " scale model with characteristic length equal to 
and N = 1, 


Some prototype trajectories and model trajectories, adjusted in the 
above manner, are compared in figure 2. The agreement is improved as 
Zq is increased. It should be noted that adjusting Zq ^ in the afore- 
mentioned manner results in slightly different trajectories and time- 
distance histories for model and prototype. Consequently, the inter- 
ference field will affect the model and prototype differently. This 
difference is believed to be negligible in most cases. 

The equations for the rotational motions are 


9 = 


C 0 § P M 2 K s L 3 g KQpg 

^22 ~~ 2 

dKyl/Kp L dL 


KePS 

0 - 0Q + £— 

dL 


K0Pg 2 

"u 


0 = 0 o t + ? 

2dL 2 


( 10 ) 

( 11 ) 

( 12 ) 


p 

The rotational accelerations vary inversely with L rather than 
with L as was the case for the linear accelerations. Substituting t x 

(from eq. (A)) for t in equation (12) gives 
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K e N 

K x 


(13) 


Therefore, ©tx is independent of L provided that L9 q and are 

constant. The same substitution into equation (ll) along with the condi- 
tion that L9q = Constant gives 


3 tx 


Constant ^9 [2 . pNg 
L L t/dKx 


(14) 


Therefore, 


tx 


is inversely proportional to 


L. 


Equation (13) establishes that 9^ x is independent of L; therefore, 
the contribution of 9^ x to the angle of attack at t = t x is correct. 
However, in model tests, z at t = t x is always too small because the 
contribution due to gravity is too small, even when the condition for 
z 0 M (eq. (9)) is imposed. Consequently, the contribution of z to the 

angle of attack at t = t x is too small. The angle of attack due to z 
at t = t x , a^. ^- x , is approximated by the ratio of z^. x to the hori- 
zontal velocity of the model relative to the air. Writing this ratio by 
means of equations (2) and (6) with the substitution of t x (from 

eq. (4)) for t gives 


a z,tx 




2gdNL 2 

V K X P 


U - 


2K x pgN 

d 


(15) 


Writing equation (15) for model and prototype conditions, subtracting, 
and inserting the condition for z 0jM (from eq. (9)) gives the following 

for the error in the angle of attack at t = t x : 



( a z,tx) M ‘ ( a z,tx) p ^x 


(16) 


L 


NACA TN 3907 


9 


As an example, using the same model values as for the example 
following equation ( 9 ), and setting U equal to 2,000 feet per second 
gives for equation (l 6) 


If the condition Zq,M = ^ 0 ,P is used instead of equation (9), the error 

will be just twice as large as that given by equation (l6) . In model 
tests where this method of simulation is employed, the angle of attack 
at t = t x is always too small by an amount given by equation (l 6). 

Should particular conditions warrant, M could be increased by an 
amount equal to Ao^. 

The period of the short -period longitudinal oscillation, P, is well 
approximated by 



'tx 


2000 


V/(2)(32.2)(1)(10) 


= -0.006l radian or -0.35° 


P = 2 n 




(17) 


The period is proportional to L provided that — is constant 

P 


The time to damp the short-period oscillation may be written as 



0.693 



(Constant) £\/TL 
P 


(18) 
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The time to damp is proportional to L provided that — , and T 

are constant. 

The conditions and results for the foregoing light -model simulation 
are summarized in table I. 

The condition that is constant results in the following equa- 

d 

tions for the model and prototype. 


"M 



(19) 

Pp\Lp/ 


dpH 

(20) 

P W 



Modified Light-Model Method 

As pointed out earlier, the gravity component of the vertical accel- 
eration is deficient in model tests utilizing the light-model simulation. 
One method of circumventing this problem is to accelerate the parent 

model upward at a value equal to - ljg« A photograph of a test setup 

utilizing this idea is shown in figure 3- The parent model is mounted so 
as to have vertical freedom only; the acceleration in this case is produced 
by a combination of wing lift and a pneumatic cylinder. With this arrange- 
ment, the trajectory of the store relative to the parent model may be 
made to duplicate exactly the prototype trajectory. Equations (l) to (4) 
and (10) to (20) apply without modification; equations (5) to (9) are 
inconsequential because the trajectories are correct without adjusting 
z Q m . The modified light-model simulation is summarized in table I. 


Heavy -Model Method 

In this method, the aerodynamically produced accelerations are forced 
to be independent of L. The equations of linear horizontal motions are 


x = 


KxPg 


= K. 


ax e 


dL 


( 21 ) 
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* = *0 + ^xSt 

(22) 

x = V + \ K ax St2 

(23) 


where the value of K ax 
is held constant during 

From equation (21) 


In this method 2. q s proportional to L. The substitution for f- 
u d 

from equation (24) will now be made for the equations in the previous 

sections of this report. 

Solving equation (23) with x Q = 0 and 

t = /Ml 

From equations (22) and (23), x tx and x tx 
and L, respectively, provided that = 0. 

The equations for vertical linear motions are 

2 ■ ( 1 + ^) g * ( 1 + M 8 (26) 

i ‘ * ( L + K a z )6‘ (27) 

1 - V + l( X + K az)st 2 (28) 


x = NL gives 


(25) 


are proportional to 


is the value of the horizontal load factor which 
scaling. 


p K ax L 
d " K* 


(24) 


where the value of K ax is the value of the vertical load factor which 
is held constant during scaling. 
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Comparison of equations (21) to (23) with equations (2 6) to (28) 
shows that the vertical and horizontal motions scale identically provided 
that x 0 = z 0 = 0. Consequently, the trajectories, normalized with respect 

to L, will be identical. However, for the case where Xq = 0 and 
z 0 £ 0, this is not true and it will he necessary to adjust to 

make (-) equal for model and prototype. Substituting t x (eq. (25)) 
\x/t x 

into equation (28) and dividing by x = NL gives 



0 


' ^axS^ 


1 + *az 

Kax 


(29) 


Writing this equation for model and prototype, equating, and solving 
for z Q?M gives 


z<D,M 



(30) 


In the heavy -model simulation, z^ ^ is less than z^ pj whereas, 

in the light-model s imulat ion , the reverse is true. An example of the 
heavy -model simulation is given in figure 4. 

Placing the restriction of equation (30) on equations (27) and (28) 
and substituting t x for t shows that z^. x is proportional to \flj 

and z^_ x is proportional to L. 

The equations for rotational motions are obtained by substituting 
for from equation (24) in equations (10), (ll), and (12). 

0 = ^ 22 ? (3D 

Kx L 


• • 

0=0 


0 


^0^ax^ 

KxL 


( 32 ) 
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e = eut + 1 K - eKa ^ S + 2 


2 y, 


Substituting t x (from eq. (25)) for t gives 


e tx = e 0 


2NL + ^ 
K ax® + K * 


(35) 


(34) 


Therefore, 0 tx is independent of L provided that x/l0q is constant. 
With the substitution for t and the condition that = Constant, 

equation (32) gives 


Constant 


tx 


& 


K x V L 


(35) 


Since 0tx is independent of L, its contribution to a^ x is 
correct; however", because z^ x is proportional to ^L, the contribution 
of z^ x to Of. x is too small. Following the development for equa- 
tion (15) and using equations (22), (25), and (27) gives 


a z,tx 



(36) 


Writing for model and prototype and including condition for z^ x 
(from eq. (30)), gives 


^“tx ( a z , tx) m " ( a z,tx)p 

*0 ,«( X - n/§) + ( L ^ K az) ^ 


(37) 


u 
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If the condition z^ M = ^0 P use< ^ 

correction is given by equation (37) except 
Substituting some typical values gives 

(20) (l - \/£o) + (l + 0.5) (j ° * 5 " 

^tx = 2000 

= -0.0452 radian or -2.6° 

This error may be appreciable, and whether it may be desired 
to adjust a Q M will depend on the particular testing situation. The 

error is larger for the heavy-model simulation than for the light-model 
simulation because z^_ x is approximately proportional to yjli for the 
former and approximately independent of L for the latter. 

The period P obtained from equation (17) by imposing the condition 
of equation (24) is 

P = (Constant) JT, (38) 

The period is properly scaled. 

The time to damp obtained from equation (l8) and equation (24) is 

t^y2 = (Constant) >/T (39) 


instead of equation (30), the 


that z 


0 ,M 


1 - 


= 0 . 


The models will therefore be too lightly damped. 

The heavy -model simulation is summarized in table I. The condition 
that £=4^ gives 

d Kx 


Lp Pm 

dw = v- — d p 

H* P P 




(40) 


( 41 ) 
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CONCLUDING REMARKS 


Two methods, the light model and the heavy model, have been presented 
for dynamically scaling store models for wind-tunnel store -dropping 
studies. In each, it is assumed that the model and prototype Mach numbers 
are necessarily equal and that Reynolds number effects are negligible. 

The light -model method gives exact simulation for all the motions except 
that the gravity component of the vertical acceleration and, consequently, 
the vertical displacement between the store model and the parent airplane 
model is deficient. The effect of this deficiency on the trajectory is 
reduced as the vertical ejection velocity is increased. This deficiency 
is eliminated completely by a modification of the light-model method, 
wherein the parent model is accelerated upward at a rate which offsets 
the gravity deficiency. This modified light -model method provides complete 
simulation of the store motion for any ejection velocity. Either of the 
light -model simulation methods will simulate accurately the oscillatory 
characteristics of the store. The basic light-model method may be used 
when the ejection velocities are large; the modified method should be used 
when the ejection velocity approaches zero. 

The heavy -model method simulates exactly the trajectory of the store 
when the parent model is stationary; the need for accelerating the parent 
model as in the modified light -model method is obviated. The heavy -model 
method gives exact simulation of all the store motions except that the 
short-period longitudinal oscillation is too poorly damped; this defect 
is inconsequential because the time during which the store model is 
critically close to the parent model is generally small compared with the 
period of the longitudinal oscillation of the store. Thus, when simula- 
tion may be achieved by the heavy -model method, it is recommended. Expe- 
rience has shown, however, that it is often impossible to make the models 
sufficiently heavy and with the proper moment of inertia to satisfy the 
heavy-model method, particularly when simulating high altitude drops in 
a sea-level facility such as the preflight jet of the Langley Pilotless 
Aircraft Research Station at Wallops Island, Va. Store construction 
problems may, in fact, dictate the method of simulation. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 3, 1956. 
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APPENDIX 

TEST METHOD USED IN THE PREFLIGHT JET 


Test-Facility Characteristics 

The tests in the test -unit of the preflight jet of the Langley 
Pilotless Aircraft Research Station at Wallops Island, Va., are made in 
a free jet issuing from approximately 27 -inch- square nozzles at subsonic 
Mach numbers up to approximately 0-95 and at supersonic Mach numbers 
of 1.4, 1.6, 1.8, and 2.0. The static pressure of the jet is adjusted 
to equal the static pressure of the ambient still air and is therefore 
approximately equal to 2,ll6 pounds per square foot. The stagnation 
temperature may be varied between 600° and 1,060° Rankine. 

At supersonic speeds, the region upstream of the nozzle exit avail- 
able for testing is contained within a wedge defined as follows: The 

base is the nozzle exit, the sides are inclined at the appropriate Mach 
angle, the vertex edge is vertical, and the ends are the flat top and 
bottom walls of the nozzle. (The nozzles are two-dimensional.) The 
region downstream of the nozzle exit available for testing is defined by 
a pyramid the base of which is the nozzle exit and the sides of which are 
inclined at the appropriate Mach angle. (See fig. 5-) Because the high- 
speed photographs are taken from the side of the jet along a line at right 
angles to the nozzle center line, only the region downstream of the 
nozzle exit can be used to observe the dynamic behavior of a model; the 
region inside the nozzle and upstream of the exit may, however, include 
those parts of the parent aircraft required for proper simulation of the 
interference flow field. 


Photographic Method 

The high-speed photographs are made by continuously illuminating 
the model and exposing the film through narrow radial slits in a 
rapidly rotating disk placed in front of a camera. Up to 1,000 exposures 
per second may be taken. The result is a series of images of the model 
taken at fixed time intervals superimposed on a single 8- by 10-inch 
photograph. The test setup is shown in figure 3; a typical drop photo- 
graph is shown in figure 6. 
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TABLE I.- SUMMARY OF SIMULATION METHODS 


H 

00 


PrhP Mach numbers air properties, static air temperatures , nondimens ional radii of gyration in pitch 
[The Mach numbers ^ air prop rp , f mode! and prototype. The models must, of course. 


and the aerodynamic coefficients must be ecmal for model and prototype, 
be geometrically similar to the prototypes .J 
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Origin of axis system at point of store 
release-, stays with parent model in 
modified light-model method. 



Figure 1.- Coordinate system. 
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(a) Space trajectories. 

Figure 2.- Comparison of model and prototype trajectories for light-model 
method. Curves calculated assuming uniform flow field; z 0 ,M ad j usted 
so that (z/x)^ = (z/x) p when x/l = 1.0; K^p/d = 10; K^p/d = 5; 

Lp = 10 feet; Im = 0.5 feet. 
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(b) Time histories. 


Figure 2.- Concluded. 



l-94715.1 

Figure 3.- Typical free-drop setup in the preflight jet. 
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(a) Space trajectories. 

Figure 4.- Comparison of model and prototype trajectories for heavy-model 
method. Curves calculated assuming uniform flow field; K ax = 1.0; 

^az = 0*5; Lp = 10 feet; 1^ = 0.5 feet; curves are coincident for 
values of i Q specified. 
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(b) Time histories. 
Figure 4.- Concluded. 
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Side View 


Figure 5 .- Region of preflight jet available for drop test. 
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Figure 6.- Typical photographic history of simulated, bomb drop. 


NACA TN 3907 



